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1 Introduction 


The fractal structure of space-time is of primary interest in any theory of quantum gravity. 
Although the starting point usually is manifolds of a specific dimension d, the quantum 
theory instructs us to perform the average over all equivalence classes of metrics and in this 
way the “effective” dimension of space-time can be different from d. Presently we do not 
have a well-defined theory of qnantnm gravity in fonr dimensions. Two-dimensional gravity 
is an interesting laboratory, where we have available both analytical and numerical tools 
with which we can stndy the fractal strnctnre, and in the following we will confine onrselves 
to two-dimensional quantum gravity, although some of the definitions and relations are 
given for arbitrary d. 

We assume that the partition function for Euclidean quantum gravity can be written 
as: 

Z{k) = (1) 

where the integration is over equivalence classes of metrics [gab]-, and the action is the 
gravitational action, nsnally taken to be the Einstein-Hilbert action which depends on the 
gravitational constant G. (j) symbolizes the matter fields, and the dependence of the matter 
conpling is snppressed. It is sometimes convenient to consider the partition fnnction for a 
fixed volume V. Since the volume is conjugate to the cosmological constant, we can define 
the finite volnme partition fnnction withont reference to A: 

Z{V) = 11 VlgabW 5{J^-V) (2) 

The partition fnnctions in (1) and (2) are related as 

fOO 

Z(A) = / dVe-^^Z{V). (3) 

Jo 

Recently it was shown [1, 2, 3] that the volume-volume correlator of Euclidean quantum 
gravity is a perfect probe of the fractal structure of space-time. It is defined as 

Ga{R) = jj d'^^^)d'^i'^f^)5{d,{i,0-R). ( 4 ) 

where dg{^,^') denotes the geodesic distance between the points labeled by ^ and 
calculated with the metric gab- One can view GA(i?) as the partition function for universes 
with two marked points separated by a geodesic distance R. If scaling argnments can be 
applied to the system, one expects a generic behavior 

Ga{R) ~ A''-^‘F{A'^R), (5) 


where jg is the string susceptibility of the system, i.e., 

ZiA) const. A^ -|- less singular terms. ( 6 ) 


The behavior (5) follows from the definitions since 


d‘^Z{A) 

d‘^A 


dRG^iR) 


~ const. A -|- less singular terms. 


( 7 ) 


However, there is a slight subtlety associated with (5) and (7). If — 1 < 7 ^ < 0, as is the 
case for nnitary theories with 0 < c < I, the coefficient in front of (7) is negative (cf. [1]). 
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Since the function F in (5) is positive for unitary theories and falls of fast for large R (see 
(18) below) this is only possible if the function F{x) is singular for small x, such that the 
term is a subdominant term, while the dominant term is analytic in A. 

The exponent 7 ^ determines the proliferation into so-called baby universes [7, 8 ] and the 
scaling exponent v can be identified with the inverse Hausdorff dimension of the system: 


V = l/dn- 


( 8 ) 


We denote (Ih the grand canonical (intrinsic) Hausdorff dimension, since it is defined in the 
ensemble of manifolds with fluctuating volume Vg = f but constant cosmological 

constant A. Eq. ( 8 ) is reasonable since the average volume of the ensemble of manifolds 
with partition function G\(R) is: 


iy,)R = 


d\nGA{R) ^ 


(9) 


for R ~ . This relation follows from scaling behavior (5). 

In the case of pure two-dimensional Euclidean quantum gravity one can calculate 
G\{R) analytically [1]: 


Ga{R) 


A3/4 


cosh 

sinh^ Ai/4i?' 


( 10 ) 


If we expand Ga{R) for small R we get: 


Gk{R) = ^ - + 0{k^R% ( 11 ) 


i.e. the first term is indeed singular for small R and analytic in A. 

We can define the two-point function Gv{R) on the ensemble of metric manifolds con¬ 
tributing to the partition function ( 2 ) in the same way as we defined Gl{R) corresponding 
to the ensemble (1). It can be viewed as the partition function for universes with a fixed 
volume V and two marked points separated by a geodesic distance R. Gv{R) is related to 
Ga{R) by a Laplace transformation, in the same way as Z{K) is related to Z{V) by (3): 


fOO 

Ga{R) = / dV Gv{R), 

Jo 


( 12 ) 


and Gv{R) is by definition related to the average volume S{R)dR in a spherical shell of 
thickness dR a geodesic distance R from a given point: 


{S{R))v = 


1 Gv{R) 
V Z(V) 


R 


F{—) 




(13) 


where F{x) is a function with -E(O) > 0 , which falls off for large x. Eq. (13) is another 
definition of the (intrinsic) Hausdorff dimension associated with the canonical ensemble of 
metric manifolds used in two-dimensional Euclidean quantum gravity. We denote d^ the 
canonical Hausdorff dimension. A priori there is no reason that d^ = dij, but in fact it 
follows from the definitions that if F{0) > 0 then d^ = dn- Clearly E(0) > 0 is a necessity 
for (13) to be valid as a definition of the Hausdorff dimension in the canonical ensemble. 
But from (5) it follows by inverse Laplace transformation and the assumption (13) with 
E( 0 ) > 0 that 

Gv{R) ~ for i? 0 . (14) 
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This is only compatible with (13) if = 1. In the following we will not distingnish 
the two Hausdorff dimensions, but completeness it should be mentioned that one can find 
non-unitary models where dh > dn- An example is the so called multicritical branched 
polymer model [9], where d^ = 2, while dn = m/{m — 1), and m > 3 is an integer. Prom 
the relation between GA(i?) and Gy(i?) we deduce that 

GAiR) ~ (15) 

where F is related to F by 

roo , 1 

F{x) ~ dye-y- (16) 

We have that P(0) > 0 since the integral in (16) is well defined for a; = 0, as F{y) is 
expected to fall off faster than a power of y (see (19) below) for y —)■ oo, and 7 < 1. In 
fact we can say somewhat more about Ga{R) since the integral after R should be as in 
eq. (7). This implies that we should be able to split Ga{R) as follows: 

Ga{R) = (17) 

where P{x) is a polynomial of degree less than [1 — 7 ^] with P(0) > 0 , while F{x) is a 
non-analytical function, also with P’(O) > 0. P will not contribute to the singular term 
in (7) and the entire contribution comes from the integral over the last term in (17). 
Under an inverse Laplace transformation the first term in eq. (17) will not contribute to 
Gy(P), U > 0, and in this way the small R dependence of the second term 

in (17) is the analogy of the term Gv{R)- It is seen that this scenario is 

indeed satisfied in for the exact solution (lO)-(ll) of pure gravity. Further aspects of the 
expansion of Gy(P) is discussed in the appendix. 

Prom general arguments we finally expect that Ga{R) falls off exponentially [4] 

Ga{R) ~ R > (18) 

and this translates by inverse Laplace transformation to 

Gv{R) ~ R > UT (19) 

While (18) and (19) are exact in the case of pure gravity, they are unfortunately only 
bounds in the case where we have matter coupled to gravity. 

A convenient quantity from a numerical point of view is the normalized distribution 

ny(P) = 1 {S{R))y = U-V'‘-ip(a:), x = ^, (20) 

such that 

poo poo 

1 = / dRnv{R) = / dxx^^~^F{x). (21) 

Jo Jo 

Eq. (20) has the form of a finite size scaling relation in the sense that it is a universal 
function of the reduced variable x times a scale factor V~'^. In numerical simulations eqs. 
( 20 ) and ( 21 ) turn out to be quite useful. 

The Hausdorff dimension is not the only dimension which can naturally be defined on 
an ensemble of manifolds. The so-called spectral dimension is defined as follows: For a 
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given manifold the propagation of a massless scalar particle is described by the inverse 
Laplacian, where 

A, = (22) 

v9 

The inverse Laplacian has a heat kernel representation, 

fOO . . 

= I dT (23) 

where the heat kernel ,T) = is the kernel of the diffusion equation 

= Aa*. (24) 

and the normalization is / Kg{^,^',T) = 1. In this paper we consider the initial 

condition, 

= {e\0 = 

To the diffusion equation we can in the usual way associate a random walk interpreta¬ 
tion. In this way K(^,^',T) is the probability per unit volume that a random walk, which 
starts at ^ will be at at “time” T. 

The heat kernel K{^,^',T) has the asymptotic short distance expansion (the so-called 
Hadamard-De Witt-Minakshisundaram expansion, see for instance [10] for a review) 

-d2(^,^')/4T 

K,{(,(',T)= AK.i'.T), (26) 

where dg{^,^') again denotes the geodesic distance and 


A(^,^',T) ~ A(^,^,0) = 1. 

?'=0 


(27) 


The functions satisfy certain invariant differential equations, and in the coinci¬ 
dence limit ^ one can express entirely in terms of local invariants, of which 

the first ones are: 


« o (^,0 = 1 , 
ai (^,0 = 

«2(^, 0 = - K67^“') + ^A,7^, (28) 

where IZabcd^ TZ are curvature tensors. By taking the trace of the operator Kg 

corresponding to the kernel Kg{(^,(^',T), we get 

Tr Kg{T) ^ jd^^^)Kg{^,^,T) ~ ^ f; (29) 

r=0 


where 

Ar = ld^^^)ar{^,0. 


(30) 
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We call the power dg the spectral dimension. For a smooth manifold we have dg = d. 
However, it is possible to define diffusion and consequently the spectral dimension dg on 
more general structures than manifolds. We will discuss such cases later. If V is the 
volume of the manifold it is seen that 

— Ti- Kg{T) = average return probability of a random walk at “time” T. (31) 

Since it is known that the kernel can be expressed entirely in terms of 

reparametrization invariant quantities, as indicated in the asymptotic expansion (26), it 
makes sense to talk about the average in quantum gravity. For the asymptotic expansion 
(29) we get 

^ 00 

(TiK,{T)) ~ (32) 

?'=0 

where the quantum gravity average is taken over the operators which enter in (28)-(30). 
Strictly speaking we cannot be entirely sure that the quantum average can be taken term 
by term and we define the quantum spectral dimension dg by 

(•& ii-,(T)) ~ ^ tor T^O. (33) 

and in case we consider the average over manifolds with fixed volumes the quantum return 
probability will be given by 

^('»Ka(0),,. (34) 

More generally we expect the behavior of (26) to be replaced by one which involves 
the geodesic distance as defined for the volume-volume correlator: 

jydh — l 

KviR,T) (35) 

where H is a fixed volume of the manifold, H{0) >0 and^ 

Ky(R,T) ^ 

X j (36) 
With this definition it follows that 

fdRKHR,T) = l. = (37) 

and that the average geodesic distance travel by diffusion over time T is given by: 

{R{T))y = / dRRKv{R,T) (38) 

■Jo 

It is easy to understand that dh can be different from d in quantum gravity. We defined 
dh by 

^im{5(i?))^ ~ R^'^-\ (39) 

®An alternative definition would be one where we divide with Gv{R) in (36), rather than with VZ{V). 
as in (36). By (13) it corresponds to the removal of the factor from Kv{R,T). 
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where the average is performed over all metric manifolds of fixed topology and volnme 
V. If we could take the limit i? —)■ 0 for each (smooth) manifold, before taking the 
quantum average, we would get dh = d. However, there is no reason such an interchange 
of limits should be possible, and as mentioned above, it is not possible for pure gravity 
in d = 2, where dh = 4. For the spectral dimension the situation is less clear. It seems 
quite reasonable that the asymptotic expansion (32) is valid, in which case the spectral 
dimension is dg = d. However, even in this case we will find some situations where the 
spectral dimension dg is different from d. 

The purpose of this article is to report on extensive numerical simulations where we 
determine the Hausdorff dimension as well as the quantum spectral dimension of quantum 
gravity coupled to various matter fields^. 

The rest of this article is organized as follows. In sec. 2 we outline the numerical setup 
and in sec. 3 we test the scaling prediction (20) for pure gravity, while sec. 4 contains 
the generalization to the theory of gravity coupled to matter fields. Sec. 5 deals with the 
analysis of the so-called spectral dimension. In sec. 6 we analyse the use of finite size scaling 
relations in quantum gravity and make predictions about the behavior of reparametrization 
invariant two-point functions. Finally, sec. 7 contains a critical discussion of the results 
obtained so far. 

2 The numerical method 

The numerical setup of the simulations presented in this paper is as follows: For all stud¬ 
ied systems we use the discretization commonly known under the name of the “dynamical 
triangulation”. We consider ensembles of surfaces built of equilateral triangles with spher¬ 
ical topology and a fixed number of triangles. These surfaces can be viewed as dual to 
planar (f)^ diagrams. In the language of the (f)^ theory we include diagrams containing 
tadpole and self-energy subdiagrams. In the language of the direct lattice it means that 
we allow “triangulations” where two vertices are joined by an arbitrary number of links 
and where a link can join a vertex to itself. The coordination number of a point can be any 
positive integer. All these diagrams are present in the simplest matrix model formulation, 
without the necessity to renormalize the couplings. In earlier works it has been noted that 
inclusion of these subdiagrams actually speeds up the convergence when critical indices 
are measured (cf. [11]). A similar formulation was used in [5]. In case matter fields were 
used in the simulations their interactions were such that they could be viewed as placed 
in the centers of triangles (or on vertices of the dual (f)^ graphs). 

A new element in simulation is the use of so called generalized minbu surgery^ in the 
update scheme of the surface geometry. These are global rearrangments of the triangula¬ 
tion which are used in addition to the the standard “flips”. They are closely related to the 
“minbu surgery moves” described in [12] for the two-dimensional system without tadpole 
and self-energy subdiagrams or in [2] for the four-dimensional systems. For the sake of 
the present simulation we define a “minbu” to be a smaller part of the triangulation sepa¬ 
rated from the remaining part by two links, joining the same two vertices (in the dual (f)^ 
language this corresponds to a special form of the self-energy subdiagram, special in the 
sense that we exclude the case when two vertices have the same label). The two links form 

■^Closely related work has recently appeared in [5]. Where overlapping, the conclusions are indentical. 
However, as we will show, the situation in pure gravity is actually much better than it appears from the 
“raw” data presented in [5]. 

^ Minbu is an abbreviation for “miminal neck baby universe” [7, 8]. 
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a “neck” of the minbu. The new algorithm consists of two steps: in the first the snrface 
is cut along the neck and the boundaries are closed. In effect the spherical surface splits 
into two surfaces both with spherical topology. In the second step we select at random a 
li nk on each surface joining two points with different labels. The surfaces are cut along 
these links and glued together along the boundaries, the chosen links becoming a neck of 
the new minbu. It is rather simple to work out the detailed balance condition for this 
type of move. In practice we put a lower limit on the size of the minbu to avoid moving 
numerous small minbus. 

Both the standard flips and the new global moves are organized in “sweeps”. In case 
of the flips they correspond to the number of attempted flips being equal the number of 
links. For the minbu surgery moves a sweep corresponds to the number of attempted 
moves of this type being equal to the number of minbus on the surface before the sweep 
(this number usually changes during the sweep). A global sweep was in all cases followed 
by 9 standard sweeps (the global moves alone do not satisfy the ergodicity requirement). 
The new move helps to reduce the correlation times and it’s effect is quite dramatic (cf. 
[12] and [2] for discussion of this point). It also gives for free a possibility to measure 
-ystr (cf. [7, 8]). The results presented in this paper corespond to system sizes 1000, 2000, 
4000, 8000, 16000 and 32000 triangles. For these systems the lower limits on the minbu 
size were chosen in the range from 10 to 20 triangles. Since the update of the system is 
“cheap” in computer time, the typical measurements were performed every 200 sweeps, 
safely above the longest autocorrelation times observed and the number of measurements 
in a typical experiment ranged between 1000 and 5000. In the analysis of the diffusion 
equation only triangulation with 4000 and 16000 triangles were used because of the large 
measurement times. In this case the measurements form the real time barrier, taking up 
to 95% of the computer time for the larger systems and we were forced to reduce the 
number of measurements. A typical experiment in this case was an analysis based on 100 
configurations, separated by 1000 sweeps. 

3 Pure gravity 

When we use the formalism of dynamical triangulations as a regularization of the theory 
of quantum gravity the discretized volume is identified with the number of triangles N: 

V = Na‘^, (40) 

where is the area of each triangle. In addition we identify the volume elements 
with the area of the triangles and the geodesic distance r between two triangles as the 
shortest path along neighboring triangles. With this definition the discretized geodesic 
distance is always an integer. This definition is only one among many possible. We are 
later going to test the results for other definitions of the geodesic distance. In the scaling 
limit they should all be equivalent. If this is not the case, it is difficult to have any 
confidence in the “continuum” results extracted from the discretized theory. 

Pure two-dimensional gravity is a good test case for numerical simulations since we 
know the exact formula for Ga{R) and consequently for In the appendix we 

outline how to find x‘^'^~^F{x), which turned out to be sligtly non-trivial. 

We denote the discretized version of nv{R) by nN{r). It satisfies 

N 

^nAr(r) = 1. (41) 

?'=0 
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It is easy to obtain n]\f{r) directly from the measurements and we can check if the data 
are compatible with the form (20), i.e. if it is possible to find a number dh such that 

N^/^'^nN{r) = x^'^-^F{x), x = r/N^/^\ (42) 


for various sizes N of the system. In principle there is an additional constant of propor¬ 
tionality a between x and r, if x refers to the continuum variable in (20). This constant 
depends on the regularization (e.g. which class of triangulations one uses in the discretiza¬ 
tion, the definition of geodesic distance etc). In the appendix we give this constant for the 
class of triangulations used in our computer simulations and our definition of geodesic dis¬ 
tance. In (42) and the rest of the article this constant has been absorbed in a redefinition 
of x in order not to make the notation unnecessary cumbersome. 

While it is possible to find a reasonable value of this way, a much better result is 
obtained by performing a shift in the values of r and N before applying the scaling relation 
(42). This is reasonable from the point of view that the shortest distances and volumes 
are lattice artifacts where we can not expect agreement with the continuum formulas. In 
this way we are led to a “phenomenological” scaling variable x 


r + a 
+ 6 ' 


(43) 


We can “derive” (43) from the data in the following way: fix d/j = 4 and fix N and 
determine for a given value of r the value of x such that the Ihs of (42) (the measured 
function) agrees with the rhs (the calculated continuum function). In fig. I we have shown 
r as function of x for different volumes N. We observe a perfect linear relation except for 
the smallest values of r and the constant a is approximately independent of N. The slope 
changes with N and we can fit it well to + b). This is shown in fig. 2, where the 

result of a common fit like (43) to the data involving N = 1000, 2000, 4000, 8000, 16000 
and 32000 triangles are shown with the factor N^/'^ + b divided out. All data are contained 
in this graph and we see that the phenomenological scaling (43) is well satisfied. 

One important lesson from fig. 1 is that a relation like (43) is not valid for the smallest 
values of r and one should simply discard these small values. 

In fig. 3 we have shown the data and the theoretical curve for d/j = 4 and an optimal 
choice of a and b. The agreement is almost perfect and the conclusion is that we already 
see continuum physics for systems as small as 1000 triangles in the case of pure gravity 
if we include simple finite size corrections like (43) A different graphical representation of 
the data, well suited for the small x region, is obtained by plotting d\ogn]\f{r{x))/dlogx 
versus x. It has the virtue that it tests clearly the consistency of dn = d^, since the value 
dff which provides us with overlapping graphs should agree with the asymptotic value of 
the curve for small x, which determines dh- We have illustrated this in fig. 4, where the 
theoretical curve is also displayed. 

As mentioned above there are many different ways to define the concept of geodesic 
distance on the ensemble of piecewise linear manifolds. One alternative definition is as the 
shortest “link” distance between two vertices. To each vertex v we assign an area element 


dAy - ^TlyCl 


(44) 


where Uy is the order of the vertex. dAy is assumed to replace the continuum d^^y^ when 
we form the average. For a given triangulation the “link-distance” and the “triangle- 
distance” used above (the link distance on the dual (;6^-graph) can differ vastly. However, 
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as ensemble average we expect that they agree up to some constant of proportionality. 
This was proven analytically for pure two-dimensional gravity in [1], Numerically it is 
also seen very clearly. 

In the next section we will study the Hausdorff dimension when matter is coupled to 
gravity. In these cases we do not know the Hausdorff dimension. Following the philosophy 
outlined above one should now determine the constants a, b and such that the distri¬ 
butions n]\f{r) can be mapped into each other as a function of x. The best dh determined 
this way from (43) would then be the candidate for a Hausdorff dimension. If we do this, 
we get 

4 = 3.9 ±0.2. (45) 

We view this result as the typical accuracy one can expect. For choices of 4 in this 
interval we get good overlaps of the distributions with acceptable x^. In order to get a 
more precise determination of 4 one would have to improve the statistics of the large-r 
tail of the distribution n(r), which in principle is possible, but in practice is very computer 
intensive since the probability of creating these elongated universes is exponentially small. 

One can use simpler charactistics of the distribution n]\f{r), like the fact that the peak 
of the distributions should scale as to determine 4- It gives results compara¬ 

ble to (45). However, we regard the procedure outlined above as more convicing, since 
the requirement that one can map the various distributions onto each other for different 
volumes is a much stronger test of scaling. 

4 Matter fields 

We can now perform the analysis outlined above in the case of matter coupled to gravity. 
We have performed extensive computer simulations for Ising spins coupled to gravity, 
three-state Potts model coupled to gravity and one to five Gaussian fields coupled to 
gravity. In the critical point the Ising model describes a c = 1/2 conformal field theory, 
the three-states Potts model a c = 4/5 conformal field theory®, while the Gaussian fields 
automatically are critical with a central charge equal to the number of Gaussian fields. 

For all these theories we can measure nvr(r) and try to determine a possible 4- For 
the theories considered so far, it seems that we have scaling according to (42) and (43). 
We have illustrated this for c = 1/2, 1, 2 and c=5 in fig.5. Each of the graphs contain 
the scaled data for system sizes 1000, 2000, 4000, 8000, 16000 and 32000 triangles. The 
results for other c’s (c = 4/5, 3 and 4) are similar. We have shown both the distributions 
n{r{x)), as well as the derivative distributions dlog n{r(x))/dlog x which, as mentioned 
above, have the advantage of displaying the small x region more clearly. 

In these cases the Hausdorff dimensions have been chosen to be 4 for c < 1, 3 for c = 2 
and 2 for c = 5, respectively, and we see as good scaling as in the case of pure gravity. As 
already remarked the internal consistency of the fitting to a single Hausdorff dimension 
(i.e dff = dh) requires that the logarithmic derivatives should converge to the chosen dh 
for small x. As is seen on the the plots on the left hand side of fig. 5 the graphs for c = 1/2 
and c = 1 indicate a slightly lower 4 (around 3.8). The same value was found to give the 
best scaling in for pure gravity, where it is known that 4 = 4.0. Also for the other values 
of c < 1 we find the best all over scaling for 4 ~ 3.8. Finally it should be mentioned that 
for c < 1 the best values of the constants a and b in the phenomenological formula (43) 

®The value of the critical point of the three-states Potts model coupled to gravity has recently been 
calculated [13]. 
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are very close to the pure gravity values. Clearly wee have some interval of where the 
fits appear to be acceptable and in this way we extract a tentative Hausdorff dimension 
for the various theories. The result of the best dh as well as errorbars of a somewhat 
subjective nature is shown in fig. 6. Prom the figure and the fact that d = 4 gives very 
good fits for c < 1 and d^ = 2 very good fits for c = 4 and 5, it is tempting to make the 
following 

Conjecture: = 4 for c < 1. d^ = 2 for large c. 

This conjecture is in agreement with old numerical results from c = —2 systems where a 
measurement of the Hausdorff dimension gave results very close to d/j = 4. It is also in 
accordance with the recent numerical simulations reported in [5] for c < 1. 

It may be surprising (but see sec. 7) that d/j = 4 for all c < 1, since there are indications 
of different scaling in pure gravity and gravity coupled to (m, m + 1) conformal matter. 
Using a Hamiltonian formulation of string field theory for non-critical strings it has been 
argued that the scaling should be governed by the dimensionless variable where 

T is the so-called proper “time” [6]. This should be compared to the known result for pure 
gravity as presented in eq. (18). If we could identify T with the geodesic distance one could 
conclude that dh = 2m for the (m, m + 1) model coupled to gravity. However, for m > 2 we 
have no such identification. It is in principle possible to have an identification T = 
in which case the conjectures from the Hamiltonian formalism would be in accordance with 
the observed dh = 4. It is natural in the Hamiltonian formalism to expect^ an exponential 
decay like (18) of the volume-volume correlator GA(i?) with respect to T. Such a decay 
and dh = 4: would lead to G\(R) ~ j g decay of G\(R) which is slower than 

exponentially, but not power like. Although unconventional, such behavior cannot be ruled 
out a priori, since it does not violate the bound (18). By inverse Laplace transformation 
we get (see also the second part of the appendix for further discussion) 

(46) 

This behavior is different from the one of pure gravity (see (19)) and can in principle be 
tested from the observed distribution n]\f{r) of the discretized theory since the prediction 
is: 

log Wat( r) ~ a;4/(2m-i) fgj. I, (47) 

where x is the scaling variable in (42)-(43). Unfortunately this difference in behavior is 
for large x and the statistics of the tail of the distribution n]\f{r) is not sufficiently good 
to distinguish between the power 4/3 for pure gravity and the power 4/5 corresponding to 
m = 3 for the Ising model. The reason is that we have to allow for unknown subleading, 
i.e. power like correction factors to (46) and that we have to consider quite large x before 
we can ignore such factors. It would be most interesting if one could prove or disprove 
(47). The statistics is much better for small x and as discussed in the appendix we expect 
the following behavior in the case of pure gravity: 

nAr(r) ~ ^1 -|- cia:^ -|- C2X^ . (48) 

While dh = 4 leads to the use of the scaling variable x as well as to the leading term in (48), 
there is no reason the next terms should agree in pure gravity and after coupling to matter 

^But we should stress that it has not yet been proven. 
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even if = 4. In the appendix we present argnments which indicate that the leading 
correction term in (48) is replaced by an term in the case of an minimal unitary 

model coupled to gravity. This power is two or less in the case of the minimal models. 
We test the sub leading corrections to n]\f{r{x)) by plotting the logarithmic derivative as 
a function of a:“: 


dlog n(r(x)) 
dlog X 


dft — 1 + acix°^ + o(a:“), 


(49) 


if n{r{x)) ~ + cia:“ + o(a:“)). The best power of a is the one where the data have 

a definite slope for x small. In fig. 7 is shows the correction term for central charge c = 0, 
c = 1/2, c = 1 and c = 5. As long as c < 1 the results agree quite well with pure gravity, 
and the power a is definitely between 3 and 5, i.e. a very small deviation, if any, from the 
result of pure gravity. In order to explain this selection of powers of the subleading term 
it is tempting to 


Conjecture: Gv{R) —t e‘^^i^sGv{R) for V ^ and 0 < c < 1. 


It would be interesting to understand the reason for this “symmetry principle”. 

Up to this point we have not presented any evidence that coupling of matter to gravity 
creates a back-reaction on the geometry associated to the metric properties as long as 
c < 1. This is in marked contrast to the situation for c > 2. Here the Hausdorff dimension 
did change (cf. fig. 6) and fig. 7 shows that the subleading small x corrections are different 
as well. The exponent a « 2 for c = 4 and 5. However, even for c < 1 the distributions, 
scaled and shifted according to (42) and (43), are not identical, as shown in fig. 8 . Note 
that the distributions agree very well for small x, in accordance with the hypothesis that 
d/i = 4 and that even the subleading exponents a in (49) are identical. Unfortunately 
we have not found a convincing parametrization of the difference between the various 
curves for c < 1 which is actually observed for x > 3. Only for large x can we use a 
parametrization like (47) due to the problems with unknown subleading corrections, as 
mentioned above. 

As mentioned in sec. 2 the string susceptibility is an aspect of the fractal structure of 
quantum gravity we get for free in our computer simulations if we use the “minbu surgery” 
update algorithm. In contrast to the Hausdorff dimension the string susceptibility 7 ^ 
shows a clear dependence on the conformal matter coupled to gravity even for c < 1 . 
In fig. 9 we have shown a measurement of jg for the various matter theories. It should 
be compared to fig. 6 . It is seen that 7 ^ —)■ 1/2 relatively fast above c = 1. All this 
corroborates on the idea that the scaling limit for c large is that of branched polymers. 
They have Hausdorff dimension 2 and 7 = 1/2. Several remarks are in order. While 
the measured 7 ^ for c < 1 agree well with the theoretical results, 7 ^ comes out too small 
for c = 1. This is a well known effect (cf. [14]) and is due to logarithmic corrections 
which are large. If they are included one gets 7 ^ ~ 0. These corrections have not been 
included in fig. 42 where we have preferred to treat all data set identical, i.e. in this case 
without logarithmic correction, which are not present for c < 1 and probably not for c > 1 
either. The results for c > 1 are somewhat larger than the previously reported results 
[14]. However, the present measurements are performed with an ensemble of manifolds 
where the triangulations (represented as by their dual (f)^ graphs) include tadpoles and 
self energy diagrams. This may explain the difference. It is somewhat remarkable that 
the values of 7^ > 0 seem consistent which a theorem [15] which states that if 7^ > 0 and 
all manifolds (after integrating over matter fields) are counted with positive weight, then 
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7s = where m > 2. 


5 The spectral dimension 

We now turn to the measurement of the spectral dimension in two-dimensional quantum 
gravity. 

We can measure the spectral dimension as defined by (33). At the discretized level 
it is natural to consider random paths between neighboring triangles. In this case the 
discretized diffusion equation can be written: 

4>{i,t + l) = (50) 

(b) 

where (ij) is one precisely for the three triangles j with are adjacent to triangle and 
zero elsewhere. The continuum normalization of $ which corresponds to the heat kernel 
is $(^,t) = K{^,0,t), and it translates to the discretized notation as follows: 

^o(^,0) =-^o) -t (51) 

The measurement of the quantum spectral dimension as defined by (33) is performed by 
generating a number of independent triangulations by the standard Monte Carlo tech¬ 
nique, and for each triangulation solve the diffusion equation (51) and measure the return 
probability $(io,i). The triangle io should be picked with even probability since it repre¬ 
sents unit area. Prom a practical point of view it is convenient to perform an average over 
some *0 for each independent triangulation generated by the computer simulations. 

As in the case of the Hausdorff dimension it is important to test if the results are 
independent of the detailed “microscopic” definition of the geodesic distance, since it is 
by no means universal. Again a simple test of the universality is obtained by studying 
the diffusion using the links instead of the triangles as possible paths, and the shortest 
link path between two vertices as a definition of the geodesic distance. As in (44) the 
area element associated with a vertex v is proportional to the order of the vertex. The 
analogy of (51) will be 

^o(^,0) = - ^o) -t <;6o(?;,0) = — (52) 

VsiO 

and the diffusion equation (50) will be replaced by: 

4>{v,t + l) = — 4>{v',t), (53) 

{vv') 

where the summation is over the n„ neighboring vertices v' to v. The results obtained 
this way essentially agree with the results obtained from (50) and (51) but are in fact 
considerable better behaved for small t as discussed below. In the following we will show 
only the results coming from (52) and (53). 

When comparing the measured return probability with the expected two points 

are important. The constant function is a normalizable solution of the diffusion equation 
on a compact manifold, and the behavior can only be valid for sufficiently small 

times on compact manifolds. For large times (f){io,t) will just be constant. Substracting 
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the constant does not really help ns, since the spectrnm of the Laplacian is discrete on 
compact manifolds, and this just means that the next lowest eigenfunction will dominate 
for large t. We expect the behavior to be the correct one only up to times which are 
of order or the inverse of the lowest eigenvalue different from zero. The second problem 
is that the behavior is not correct for small t either, due to the discretization. 

Clearly it makes no sense to talk about t < 1, but the situation is slightly worse, as is 
illustrated by considering a discretized random walk in one dimension. In this case the 
return probability is identical zero for odd discretized times t! Here we observe for small t 
a marked asymmetry between odd and even t, and it is considerable more pronounced and 
last for longer time if we use the return probability coming from (50)-(51) rather 

than if we use (f){vo,t) from (52)-(53). 

In principle these problems should disappear in the continuum limit, since the eigen¬ 
value density increases with increasing discretized volume®However, in order to use 
as efficiently as possible the small t behavior we have imitated a continuous time in the 
following way: We can write the solution to (53) as follows 

cp{v,t) = {l + A^m,Mv',0), (54) 

where the V xV matrix {V being the number of vertices in the triangulation) is the 
discrete Laplacian corresponding to the diffusion equation (53) 

a:I7«o') = Tew»')-«<')) (“) 

^ (vv') 

where the summation again is over the vertices v' which are neighboring v. The solution 
to the continuous diffusion equation (24) is given by 

= I (e,e')^o(C',0). (56) 

We copy this formula by replacing the solution (54) by 

^{v,t) = (1 + -^A'^®‘='')))„,(;6o(w',0), (57) 

n 

where n > t. For n —)■ oo the operator in (57) obviously goes which is what we 

want. Prom a practical point of view it has the advantage that simply evolving the original 
diffusion equation (53) n steps we know the vectors [(1 + A^®‘^'')"(;fio](u), and this is all we 
need in order to calculate (l)(v, t) for t < n. This imitation of continuous time behavior 
(but with a discrete Laplacian yields a much smoother small t behavior, and the 

data in fig. 10 and fig. 11 are obtained this way. 

In fig. 10 we have shown the typical return probability as a function of time for c = 0 
(pure gravity), c = 1/2, c = 1 and c = 2, 3,4, 5. We observe that the spectral dimension is 
consistent with 2 for c < 1 and that it decreases for c > 1. 

Since we solve the diffusion equation on each of the generated manifolds we can in ad¬ 
dition determine the average length of diffusion at time t and compare with the theoretical 

®While the eigenvalue density increases at the discretized level, the continuum spectrum should be 
(approximately) constant (and of course discrete) since the continuum volume V = Na2, a being the 
lattice spacing, is assumed to be constant. It comes about because we have to multiply each discrete 
eigenvalue of the discretized Laplacian with 1/a^ in order to approximate the continuum Laplacian. 
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formula (38). The results are shown in fig. 11 for c =0, 1, 3 and 5, and they are consistent 
with a dependence 

{r)t ~ ( 58 ) 

This lead us to 

Conjecture: dg = 2 for the central charge c < 1. dh = 2ds for all finite values of c. 

6 Finite size scaling 

Finite size scaling has been a very important tool in the analysis of critical systems in 
statistical mechanics. Likewise it has been a convenient tool in the analysis of computer 
simulations of quantum gravity [16]. However, in the latter case the theoretical basis is 
not well understood. It is the purpose of this section to present the present situation in 
two-dimensional quantum gravity, where we have available some analytical tools. 

One basic assumption in the theory of critical phenomena is that of a divergent corre¬ 
lation length. Let us consider a d-dimensional spin system where the inverse temperature 
is denoted /3, the inverse critical temperature (5^ and 

t = Wc - (d) /Pc- (59) 

We denote the correlation length which diverges for f —)■ 0 by ^(f): 

^{t) ~ for f —)■ 0. (60) 

The scaling hypothesis states that ^ is the only relevant intrinsic scale. Assume that the 
singular part of the free energy per unit volume /(/3) behaves like If we are slightly 

above the critical temperature, a fluctuation away from the ordered state by ^ produces 
an increase A/ ~ The probability of such fluctuation is and it becomes 

small for A/ > 1. By this heuristic argument one arrives at ~ i.e.: 

f(t) = and a = 2 — vd. (61) 

Assuming (61) is called the hyperscaling hypothesis, a is the exponent of the singular part 
of the specific heat c(/3), since we get the specific heat by differentiating /(/3) two times 
after fd. Differentiating one time we get the singular part of the internal energy e(t), i.e. 

e{t) ~ c{t) ~ (62) 

The spin-spin fluctuations will be long ranged when f —)■ 0. Prom general arguments one 
expects the following behavior: 

~ (63) 

where g(0) = 1 and g{x) falls off exponentially for large x. The second derivative of the 
free energy per unit volume with respect to an external magnetic field is the magnetic (or 
spin) susceptibility 

X(t) = j df^x {(l){x)(l){Q)) ~ (64) 
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where (l)(x) denote the continuum limit of the spin field, i.e. we have 

X{t) ~ and 7m = u{2 - rj), (65) 

where 7 ™ is the magnetic susceptibility. These standard scaling relations can partly be 
derived by renormalization group arguments and they are readily converted into finite size 
scaling relations by assuming that the peak in specific heat and spin susceptibility occur 
at the so-called pseudo critical point where 

~ y-i, (66) 

V denoting the bulk volume of the system. It just tells us that the correlation length is 
of the order of the linear extension of the system and that singular part of the free energy 
F{t) = Vf{t) is of order one at this point. By substituting (66) in the formulas above 
we get that the singular parts of the specific heat and spin susceptibility at the pseudo 
critical point behave as: 

X{V) ~ c{V) ~ (67) 

The set of finite size scaling relations is most readily derived by assuming that certain fields 
and composite operators here the spin field and the energy operator, have well-defined 
scaling dimensions 

Af) = (d - 2 + r]i)/2 (68) 

close to the fixed point. If this is the case, we get for dimensional reasons: 

In the case of the spin system this allows us to express the critical exponents in terms of 
the scaling dimensions of the spin field and the energy density Ai°\ In the case of 
the energy density we get e{V) ~ while for the magnetization we get m(V) ~ 

/d A factor IjV has been divided out compared to (69), since the quantities e(V) 
and m{V) are defined per unit volume. For the specific heat and the spin susceptibility 
we get: 

c{V) = ^ (70) 

V Jv Jv 

x{V) = y(df^x(d%{^,{x)^,{y))^V^-‘^^'°'l^. (71) 

V Jv Jv 

If we compare with (67) we can write: 

a = vd{l-2A.f^/d), 7m = i^d(l-2A(,°V(^)- (72) 

Finally the relation of the Ae°^ and A^^^ to the renormalization group is established by 
the well-known fact that the relevant eigenvalues under a scaling with A are given by 
and where 

^ = d-A(°), y, = d-A(°). (73) 

The scaling ansatz of DDK [17] generalizes (70)-(71) to 2d quantum gravity by simply 
replacing /y d?x by f d?Xy/g^ moving the ket and bra outside the integration in order to 
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include the gravitational average over 2d-nianifolds with volume V. For a conformal field 
theory with central charge c the change in is given by: 



This is the finite size scaling according to DDK and in this way it is possible to calcu¬ 
lated the gravity modified exponents a, (5 and 7 ^ for the spin system as long as c < 1. 
What is missing in order to get a description as detailed as that of ordinary statistical 
mechanics is the concept of a divergent correlation length and a corresponding relation 
to the renormalization group. The concept of geodesic distance, as introduced in the last 
section, provides a natural framework in which it should be possible to have relations like 
(63) and (64). A possible definition of the two-point function for a given field (j) and fixed 
volume V is® 

{cl^ir)m)v = llv[g,,]Vcl^e-^ f j ^(d(^,^') ”0, (75) 

where S denotes the combined action of gravity and matter fields and the Z{V) denotes 
the partition function for finite volume V. 

Pi-om dimensional arguments we have to conjecture that the quantum gravity gener¬ 
alization of (63) and ( 68 ) should be 

{<l>ir)m)v ~ (76) 

where is the Hausdorff dimension and ^(t) the correlation length, measured in geodesic 
distance unit. For f —)■ 0 we assume that ^{t) diverges for infinite V, while /(O) is finite. 
This conjecture leads to the correct scaling 

/.V^/dh 

/ dr{4>{r)4>{t)))y^V^-‘^^^l'^ (77) 

Jo 

at the critical point. If we assume that ^(t) divergences as t~'^ we get in addition the 
analogue of Fishers scaling relation: 

Tfli = i"4(l - (78) 

where 7 ^ is the susceptibility exponent for the correlator, defined as for the ordinary 
spin-spin correlator. It should be compared to (72). d has been replaced by dh when 
combined with u and the “bare” scaling exponent by the dressed one, related as in (74) 
in the two-dimensional case. 

The conjecture (76) has the appealing feature that it attributes a physical interpreta¬ 
tion of the scaling dimension at the level of correlation functions, precisely as for a fixed 
manifold, but the question is whether it is correct. It has not yet been proven that there 

®It is possible to adapt several, slightly different definitions of the correlation functions in quantum 
gravity. The one chosen here includes an angular average. We could choose to divide by this angular 
average, either inside the functional integration in (75), or outside the the functional integration, in which 
case one (by (13)) should replace the normalization VZiV) by Gv{R). In addition one should strictly 
speaking also define the connected part of the correlator. For the scaling arguments presented here it is 
not necessary. 
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exits a divergent correlation length ^(t) in qnantnm gravity in the way defined above. 
We view this as one of the most interesting unsolved problems in two-dimensional quan¬ 
tum gravity. In principle the answer can be found by numerical simulations, although it 
seems difficult presently to measure the spin-spin correlation functions with sufficiently 
accuracyThe Ising model coupled to quantum gravity seems to be an ideal test model 
for the conjecture (76). 

7 Discussion 

In this paper we have presented arguments in favor of a spectral dimension equal two in 
quantum gravity as long as c < 1. In the same regime of central charge the Hausdorff 
dimension is measured to be close to four. These measurements are in agreement with 
recent independent numerical simulations for 0 < c < 1, as well as somewhat older sim¬ 
ulations for c = —2. In addition we have shown that both the Hausdorff dimension and 
the spectral dimension decrease for large c, while their ratio stays approximately equal to 
two. 

The observation that the Hausdorff dimension, within the numerical accuracy of these 
simulations, is the double of the spectral dimension calls for a simple geometrical expla¬ 
nation. Does the fact that we have a non-zero probability for baby universe creation per 
unit area automatically imply dh = ‘2,ds‘! 

The spectral dimension and the Hausdorff dimension seemingly test two different as¬ 
pects of the dimensionality of an ensemble of manifolds. As argued in the introduction 
it is not so easy to understand why the spectral dimension of the ensemble of manifolds 
should be different from the dimension of the underlying manifold, while we do not have 
the same conceptional difficulties with the Hausdorff dimension, but the fact that both 
seem constant for —oo < c < 1 indicates that both should be considered as intrinsic prop¬ 
erties of two-dimensional quantum gravity, independent of the coupling to matter. The 
change in spectral dimension (and in Hausdorff dimension) for c > 1 thus indicates a 
drastic change in the theory, and it corroborates on the idea that for c > 1 (or at least 
for c sufficiently large) the interaction between matter and gravity is so strong that the 
two-dimensional surface is torn apart. In the regularized, well-defined theory, described 
by dynamical triangulations, this seems to take place simply by a change in the ensemble 
of triangulated manifolds. The probability of picking a piecewise linear manifold, which is 
so degenerate (branched ?) that it does not qualify as a genuine two-dimensional manifold 
becomes one. This is seen in the computer simulations. The spectral dimension of all the 
individual triangulations decreases below two. However, it should be emphasized that the 
new theory for c > 1 seems not to be arbitrary. It keeps for instance the ratio between 
dh/ds ~ 2. For large c the theory probably degenerates to that of branched polymers, but 
there might be a region 1 < c < cq where we have a non-trivial statistical theory. It is an 
interesting theoretical problem to understand this region. 

^°In [5] the spin-spin correlation function is measured, but no power law observed. The decay is fitted 
better to an exponential decay. It is presently difficult to give an interpretation to this result. 
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8 Appendix 

8.1 Numerical calculation Gv{R)- 

The Green function Gy(i?) is the inverse Laplace transform of GA(i?), i.e. 


fC+ioD 1 / f? \ 

where c is a positive real number and U{x) = x^’^~^F{x) is the dimensionless functions 
defined by (20) and (21). The constant l/12y^ is introduced to ensure the correct nor¬ 
malization (21) of U{x)\ 


U{x) 


12^^ 


(‘C-\-ioo 


ds g 3^4 cosh(rES^/^) 
27ri sinh^(a:s^/^) 


(80) 


The discretized version ni\f(r), as defined by (41), is related to U(r/N^^^) in the limit of 
large N by 


nN{r) 


a 

A^V4 


U{ 


ar 

AG/4^’ 


(81) 


where a is a constant parameter which depends on the regularization. In the case of the 
dynamical triangulations used here (which corresponds to one-matrix model with cubic 

potential) a = ^6/(12 -|- 13v^)- 

We now calculate the function U{x) by performing numerically the inverse Laplace 
transformation. Four different methods are used: 


i) direct numerical integration after s, 

ii) analytic integration, after the mode expansion, 

iii) saddle point integration, after the mode expansion, 

iv) analytic integration after Taylor expansion around x = 0. 

Method i) works well expect for x close to 0 (i.e. in practice for x > 1.8), since a 
rescaling of s with 1/a: leads to the term e*/* in the integrand and it oscillates wildly for 
a: —)■ 0. 

Next, we consider the methods ii) and iii). The expansion 

CXD 

coshs/sinh^ s = 4 ^ n? exp(—2ns) 

n=l 

can be given an interpretation as a “mode expansion” [1], and we can write 

oo 

U{x) = 48^/n''^n‘^u{nx/2), (82) 

n=l 


where 

rc+ioo 

u{t) = / -—: exp(s — 4ts^/^) 

Jc—ioo 27r? 


( 83 ) 


Our task is now to evaluate (83) 


and next perform the summation (82). 
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It is possible to express u{t) in terms of “known” fnnctions since the inverse Laplace 
transform of exp(—(m = 1, 2, 3,...) can be expressed by the imaginary part 
of McRobert’s £^-functions. We have 


u{t) 


1 8 9^. 

^^5/2i7''L^^4’4’4’ 4 • 
3 pJ.l 1 3. 
4r(l/4) ^ 4’ 2’ 4’ 


A7rj.4 


-t^) + 


n] 

5r(i/4)t2 


2V2'. 


l-f’S 


TT 


5 5 3 7 

'4U’2’4’ 




9 3 5 3 4 

4’ 4’4’ 2’ * 


(84) 


where the hypergeometric function 1 F 3 is defined by 


1^3(0; 61,62,^3;^) 


CXD 


E 


{a)k 

{hi)k{h2)k{h)k kV 


(85) 


and (a)jfc = r(a + k)lV{a). 

In principle we can use (82),(84) and (85) to calculate U{x) for all x. Prom a practical 
point of view it works well except for the smallest value of a;, where we have to include 
very many terms in (82) and for very large x where we have to include many terms in the 
hypergeometric series (85), i.e. in practice in a the region 0.02 < a: < 10. 

For large values of x it is sufficient to calculate the integral in (83), defining u{t) by 
the saddle-point method iii). In this case we get: 


u{t) 



00 

3W3^5/3g^^^-4/c/3^ 

k=Q 


( 86 ) 


where 


Co = 1 , Cl 


115 

144’ 


C 2 


20518225 

10319560704’ ^ “ 

162818713432375 

184884258895036416 ’ 


695 47755 

“ ~ 41472’ ^ “ ”17915904’ 

1316690375 36815854075 

1486016741376’ ® ” ”1283918464548864’ 

366346837418687125 

® “ ”212986666247081951232’ '■■■ 


(87) 


In general the saddle-point series is only an asymptotic series, but if we cut the series in 
( 86 ) at A: = 5 we get very good agreement with the use of the McRobert’s £^-functions 
for X G [1.6,10]. For large x it is inconvenient to use the McRobert’s £^-functions, but 
for x > 1.8 we have good agreement agreement between the direct integration and the 
saddle-point approximation, which is of course the most convenient to use for a: —)■ 00 . 

The McRobert’s £^-functions becomes impractical for x < 0.02. For small x we can 
turn to method iv). If we Taylor expand the integrand in (80) for small x we get the 
following asymptotic expansion: 


U{x) 


-48a:^ ^ 


( 2 )^ C( ‘^k 5) 


A ;=0 


A:! 


16’ 


( 88 ) 


Eq. ( 88 ) is also obtained by substituting (85) into (82) with (84) and using the zeta- 
functional regularization, The convergence radius of the expansion is 

zero, but if we cut off the summation at k, and define this function as Uk(x), we expect to 
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get a good approximation for small x. By comparing with the results from the McRobert’s 
£^-function we get good agreement for x G [0.02,1.6] if we use Ug{x). Of course Ug{x) gives 
an excellent approximation for 0 < rr < 0.02 where we could not use easily the McRobert’s 
functions. In this way we have managed to cover the whole interval 0 < x < oo, and 
except for x G [ 0 , 0 . 02 ], by at least two independent approximations. 

8.2 Asymptotic expansions in the general case 

As discussed in the introduction the generalization of (79) and (80) is: 

rc+ioo fix D 

GviR) = / ^ Gj,iR) ~ (89) 

where 11 = l/dh and 

fc+ioo 

U{x) = const. / -—: e^s‘^^^‘F{xs‘^). (90) 

Jc—ioo 27r?. 

The leading term in the expansion of F{u) is i.e. the integrand starts with a term 

in s which does not contribute to the Laplace transformation for finite volume (it gives a 
^-function in the volume). In fact we know that the first term which can contribute must 
be the term which has to be present in the expansion. If we expand the integrand 

in (80) for small x, it contains odd powers of x. However, for half of these powers the 
function which multiplies e* is just an integer power of s, which does not contribute to the 
inverse Laplace transformation for finite area. This is the reason the expansion ( 88 ) jumps 
with powers of 4 in x. If we accept from the numerical experiments that the Hausdorff 
dimension d/j = 4 we have to use = 1/4 in (90). If we at the same time use 7 ^ = —1/m, 
for the (m, m + 1) conformal field theory, the function F{u) has to be quite special. In 
fact it is natural to expect an expansion: 

CXD 

F(u) = Y. (91) 

A;=0 

where the first k which gives a contribution to the inverse Laplace transformation for 
R > 0 is determined by the requirement that the power of u should be — 1 (in order to 
give i.e. 

— — = i.e. k = m + l. (92) 

The above scenario is not very natural, but it is nevertheless realized in the case of pure 
gravity (m = 2). For m 7 ^ 2 there are no obvious reasons to expect that the first correction 
term to U(x) should be x’^ as is the case for m = 2 since the most naive choice of correction 
term from (91) appears to be . 

There is another natural extension of pure gravity in the asymptotic form. We may 
expect 

OO 

F(^,) = ^,«^p(n)e-"“'. (93) 

n—\ 

In the case of pure gravity, a = 0, 6 = 1, and p{n) = r?. If d/j =4 for any m, i.e., 
V{x) ~ + ... for rr « 0 , we find that 

2 °° 

a = 3(1 — 6 ), h = —, and /?(n) = ^ajn^*, (94) 

^ i=l 


21 



is a natural extension. Then, the leading term of F{u) becomes 

F{u) ~ 

After the inverse Laplace transformation we have 

U{x) ~ exp(—const. 


( 95 ) 


(96) 
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Figure 1: The “shift” in r = ca; — a as described by (43) for various size systems {N =4K, 
8K, 16K and 32K) in pure gravity. The slope can be fitted to c = + 6, i.e. the 32K 

system corresponds to the graph with the steepest slope. 



X 


Figure 2: The function (r + + b), a = 5.5, b = —.45, plotted against x (i.e. eq. 

(43)) for pure gravity and for N = IK, 2K, 4K,..., 32K. 
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Figure 3: The scaled distributions for pure gravity for systems of sizes IK, 2K, 4K,...., 
32K triangles, as well as the theoretical distribution (the fully drawn line) as a function 
of the scaled variable x = (r + 5.5)/(A^^/^ — 0.45). 



Figure 4: The derivative dn{r{x))/dx plotted against x for pure garivity. The same sizes 
of systems as in fig. 3 and again the fully drawn line is the theoretical curve. 












Figure 5: The scaled distributions n{r{x)) as well as the logarithmic derivative 
dlogn{r{x))/dx for c= 1/2, 1, 3 and 5. In constructing the mapping (43) we have as¬ 
sumed d/j = 4 for c = 1/2 and c = 1, d/j = 3 for c = 3 and d/j = 2 for c = 5. The constants 
a, b determined this way are quite close to the pure gravity values for c = 1/2 and c = 1. 
In all cases the data include the following discretized volume sizes: N = liF, 2K, 4K, 
...,32K 
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Figure 6: The Hausdorff dimension as determined by the finite size scaling relation for c = 
0, 1/2, 4/5, 1, 2, 3, 4, and 5. Contrary to fig. 5 has been treated as a free parameter in 
these fits and the dots denote the best values of d^- 
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Figure 7: d\ogn{r{x))Jdx as in (49) versus rr" for c = 0, 1/2, 1 and 5. For c < 1, a « 4 
seems to give a well defined slope for a; —)■ 0. For c > 1 it seems that a « 2 if a slope at 
rr —)■ 0 shall exist. All curves include volumes N= IK, 2K, 4K,...,32K. 







Figure 8: The scaled distribution for the various conformal theories with c < 1 coupled 
to gravity. The size of the systems is = 16000. The values of c is 0 (the highest peak) 
1/2, 4/5 and 1 (the lowest peak). The parameters dfi,a and b nsed in (43) been chosen 
identical to the ones used in pure gravity {dh = 4, a = 5.5 and b = —0.45), but even if 
these parameters are chosen from the best fit to the finite size scaling formulas the graphs 
shown will almost be unchanged. Note how well the results agree for small x (a: < 3). 
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Figure 9: The measured string snsceptibility versns central charge. No logarithmic cor¬ 
rections are used in the fits, even for c = 1. This is the reason the result for c = 1 is not 
the correct one (i.e. 7 s (c = 1) = 0). 
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Figure 10: —2dlogTr K{t)ldlogt (« dg) versus i for c = 0 (top curve), c = 1/2, c = 1, 
c = 3 and c=5 (bottom curve) theories coupled to 2d quantum gravity. The size of the 
systems is A^=16K. 



Figure 11: {r{t))‘^ versus t for various theories of matter coupled to gravity (c =0 (bottom 
curve), c = 1, c = 3 and c = 5 (top curve)). The system size is = 4K. Straight lines (as 
observed) indicate 2dh/ds = 4 according to (38) 






